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Abstract 

The first terms of the general solution for an asymptotically flat sta- 
tionary axisymmetric vacuum spacetime endowed with an equatorial sym- 
metry plane are calculated from the corresponding Ernst potential up to 
seventh order in the radial pseudospherical coordinate. The metric is 
used to determine the influence of high order multipoles in the perihelion 
precession of an equatorial orbit and in the node line precession of a non- 
equatorial orbit with respect to a geodesic circle. Both results are written 
in terms of invariant quantities such as the Geroch-Hansen multipoles and 
the energy and angular momentum of the orbit. 
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I Introduction 

The construction of new exact solutions of the Einstein vacuum field equations 
describing an asymptotically flat stationary axially symmetric spacetime has 
increased in the last decade due to the several methods of generation of solutions 
from a given one, such as the Backlund transformations, the inverse scattering 
method and the HKX tranformations (cf. for instance pQ and references quoted 
therein.) Nevertheless we are very far from being able to implement the physical 
behaviour of an exact solution at will and in order to obtain results which could 
be tested experimentally we are led to use approximate expressions for the 
metrics with the desired physical requirements. 

In this paper we study the influence of the first Geroch-Hansen multipolc 
moments of order higher than three in some astrophysical situations with 

stationary and axial symmetry. As it is well known, these moments can be 
calculated from the coefficients of the power expansion of the Ernst potential 
on the symmetry axis 0], being linear the relation between both families of 
constants until the third moment, that is, the octupole. The subsequent expres- 
sions for the multipolc moments of order higher than three become considerably 
more complicated as the order increases and there is not even a closed formula 
for calculating all of them. The purpose of this paper is therefore to show to 
what extent these nonlinearities affect the motion of test particles tracing their 
orbits around a non-spherical, in principle rotating, compact mass distribution. 
Of course, these terms are irrelevant for our solar system calculations, but they 
are meaningful for highly relativistic astrophysical objects, such as pulsars, as 
stated in 0. 

With this aim in mind we calculate in section [H] the first seven terms of a 
power expansion of the Ernst potential with arbitrary values for the multipole 
moments and construct the corresponding approximate metric. This result is 
used in section IIIII to obtain an expression for the perihelion precession of an 
equatorial trajectory and in section llVl to produce the corrections to the Newto- 
nian precession of the nodes of a slightly non-equatorial orbit with reference to a 
close neighbouring geodesic circle. There is some previous work on this subject 
in [Hj and 0, but these references deal only with terms up to the quadrupole 
moment. The results will be discussed in section Ivl 

II Calculation of the metric 

The metric of a stationary axially symmetric vacuum spacetime can be written 
in a canonical form in terms of the Weyl coordinates, 

ds 2 = -f{dt - Ad(f>) 2 + j{ e ^(dp 2 + dz 2 ) + p 2 d<p 2 }, (1) 

where t and 4> are the coordinates associated with the commuting Killing vectors 
dt and d$, and the functions /, A and 7 depend only on the coordinates p and 
z. 
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The whole set of Einstein equations can be shown jS] to be equivalent to the 
following system of partial differential equations, 

s = f + *X (2) 



1 2 

T - PP + -e P + e zz = — — (e P 2 +£ z 2 ) (3) 
p e + e 

Ap 

(e + e) 



A ^ = -jr^x P (5) 

7 " = [e + e) 2 < " epEp ~ ® 
(e p e z + e z e p ). (7) 



' ( £ + e) 2 

It can be shown that the integrability of the last four equations is guaranteed 
if equation J3J|, the Ernst equation, is satisfied. Therefore, in order to obtain a 
solution of the Einstein equations with two Killing vectors, it suffices to solve 
the Ernst equation and then calculate the metric functions by quadratures. 

It is also usual to write the Ernst equation in terms of another potential, £, 
related to the previous one by the following relation, 

which satisfies another partial differential equation, 

(£f - 1) (£p P + U P + £,.) = H (£„ a + C, 2 ), (9) 

which is also satisfied by This latter form is the one introduced originally 
in [S] and allows a simple integration of the Kerr metric. 

This form is particularly useful for calculating the multipole moments and 
it will be the one we shall employ. 

In order to calculate the approximate solution, we shall write the Ernst 
potential £ as a function of two coordinates r and 8 related to the Weyl ones by 
the following transformation, 

p = r sin 6* z = r cos8. (10) 

We can implement the requirement of asymptotic flatness by writing £ as a 
formal inverse power expansion in the pseudospherical radial coordinate r, 



£ = X> r ~ n ' ( n ) 
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where the functions depend only on the coordinate 6. 

Since we are interested in having a solution which is symmetric with respect 
to the equatorial plane 9 = it/2, we shall require that the functions of odd 
order be real whereas those of even order will be taken to be imaginary. 

With this information at hand we can now proceed to calculate the metric 
functions. The function / is just the real part of e and from our knowledge of £ 
we can calculate the first eight terms of its expansion in the radial coordinate, 



n=l 

A = -2 mo (13) 
f 2 = 2m 2 (14) 



( 4 ™0 3 q 2/) /,,-n 

.13 = m 2 om2 cos (15) 



2 4 

fi = —2 to to 2 H h (6 to to 2 + 2 toi 2 ) cos 2 8 (16) 



3m 4 „ 2 8 too mi 2 4to 5 
+ 2to ™2 + — — + 

+ ( - 6to 2 to 2 - Mz^giT) cos2 fl - 35 ^; os4 g ( i7) 



3TO TO 4 ^m 2 2 4m 2 TO 3 I6m 2 mi 2 ^<im 



35 45 

r. 

~35 

/35to to 4 9to 2 2 \ 4 
+ h 10 Tti\ 7713 H — cos 4 6> (18) 



f „ 2 s 356too 2 to,i 2 \ 0/1 
+ —15 too to 4 — 6 toi m 3 — 3 to 2 + 4 to 2 to ' H — cos # + 



5to 6 3m 2 m4 4to toito 3 2 16mi 2 m 2 

/7 = "8 2 U ^ "^31— + 

2TO 4 m 2 32too 3 TOi 2 8to 7 
+ 3 + 63 315~ + 
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i iOSme 2 216m TOiTO 3 2 

+ | 3 h 15 m m 4 H — h o m m 2 - 

8 11 

38rai 2 m2 n 4 1208mo 3 mi i 



— 2 m n m 2 cos 6> + 

11 02 105 J 

315 35mo 2 ?ri4 340 m mi m 3 114m 1 2 rn 2 

8 2 ll ll 

- 9 m m 2 2 ) cos 4 0-^^ cos 6 0, (19) 

8 

where the constants m n which arise from the integration of the Ernst equation 
are real if n is even and otherwise imaginary. 

The first six terms of the metric function A can be obtained by direct inte- 
gration of the equations l@J and J5j. In spite of the factor i before the expression 
for A, this function is obviously real since the constants m n are imaginary for 
odd n, 



A = -isin 2 eJ2 A nr' n (20) 

A l = -2m 1 (21) 
A 2 = -2m a m 1 (22) 



8m 2 m 1 2 
A 3 = m 3 5 ?7i3 cos (23) 



3moW3 mi m 2 16 mo 3 to 



15 



+ 



15to.qto,3 3toiTO 2 



cos 2 6 (24) 



A 



5 



3«5 4mo 2 TO3 8mo?7iiTO2 64mo 4 TOi 

+ 3 + 7 105 
21 m5 20mo 2 TO,3\ 2 63m,5 4 



cos 2 - rip cos 4 (25) 
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40mo 3 m3 8 mo 2 m\ 7772 16momi 
m 2 m 3 1 — | cos" 

105 mo 7775 35 mi 77i4 hvnim^ 



.,2 , 



cos 0. (26) 

The only function which remains to be calculated is 7 and can be obtained 
as a quadrature from equations JBJ and Q, 



7 = E^ 2n ( 2? ) 



71 = -22^ (28) 



72 = 



3 mo 777,2 mi / 9 mo 777,2 5 mi \ 2 



4 4 \ 2 

15 7770 777,2 9 mi 2 



cos 4 (29) 
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T3 = 1 ■ 1 h 
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—8 2~ + ~ 35— )C0S " 

175 777o m4 55 mi 7773 117 7712 2 . . 

cos 



74 



8 2 8 

105mom4 35 mi m3 75 m2 2 
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35momg 15mim5 45m2m4 
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9m 3 2 14mo 2 mim 3 2momi 2 m 2 16mo 4 mi 2 
~~32 165 33 + 1575 

245mom6 135mim5 315 7772 7774 81 7713 2 



cos 2 



cos b (30) 
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+ 



1617 rrio m e 1071miTO5 2415to-2TO4 625 



16 
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8 

2079 mi TI5 
32 



16 
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64 



8 

1225 m 3 2 



cos 6 9 



32 



cos 8 9. (31) 



III Perihelion precession of a closed orbit 

It is well known from Bertrand's theorem (cf. for instance [§]) that stable 
bounded orbits of particles moving under the influence of a central force which is 
neither Newtonian nor harmonic are not closed. Therefore whenever the source 
of the gravitational field is not exactly monopolar, the bounded trajectories on 
the equatorial plane will no longer be the elliptic orbits described by Kepler's 
first law but will take the form of a precessing ellipse if the deviation from 
spherical symmetry is small. 

The situation becomes a bit more complicated in general relativity. Al- 
though it has been proven ^H] that there are just two asymptotically flat static 
spherically symmetric spacetimes in which the stable orbits are closed, they are 
rather different from the classical physical situations. Therefore, when rela- 
tivistic effects are taken into account, not even the motion around a spherical 
distribution of mass is closed. This effect has been tested in our solar system 
and it amounts to a slow precession of the perihelion of the orbit of Mercury. 
Of course other multipole moments of the mass distribution will also contribute 
to this effect and, in principle, these moments could be calculated by measur- 
ing the precession of a certain number of test particles orbiting at conveniently 
different distances from the gravitational source. 

If the test particles are small enough for the tidal forces to be unimportant 
within the characteristic length of the particle, we can regard them as point 
particles. If the effects due to their intrinsic angular momentum can be taken 
as negligible it can be assumed that they trace out timclike geodesies in the 
spacetime surrounding the gravitational source. Hence, in order to study the 
influence of the far field multipole moments of the gravitational field, we shall 
have to solve the geodesic equations for the previously calculated metric. We 
shall restrict ourselves to the equatorial plane 9 = tt/2. 

Since the timelike and azimuthal coordinates are ignorable, we have two 
first integrals for the motion corresponding to the conserved quantities E and I, 
respectively the total energy per unit of mass and the projection of the angular 
momentum on the z axis, also per unit of mass, of the test particle. In terms of 
its 4- velocity u — (t, f, 9, 0) these quantities have the following form, 

E = -d t -u = f(i- A4>) (32) 

l = d^-u = fA(i-A<j)) + jr 2 <j), (33) 

where the overhead dot stands for the derivative with respect to proper time. 
Therefore the equations for t and (f> can be written as follows, 
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■ E r J- EA 

t= J + fA—^. (35) 

Another integral arises from the fact that the trajectory is timelike and 
therefore u ■ u = — 1. For the geodesies under consideration this means, 

- 1 = -/(* - A<\>) 2 + j(e 2 -' r 2 + r 2 ft). (36) 

From the previous three equations r can be obtained as a function of the 
non-ignorable coordinates and the conserved quantities. However, since we are 
interested in the shape of the orbit rather than in its time evolution, we divide 
<|36l) by <f> to get the derivative of the radial coordinate with respect to the 
azimuthal angle, 

r, 2 -e- 2 A^— ?L-A (37) 

It will be useful to write this equation in terms of another function u = 1/r 
as it is done in classical mechanics for solving the motion under central forces, 



V = e ~ 27 ( /2( f_/ A)2 - « 2 i = = E c « u " + °(" 7 )- ( 38 ) 



n=0 



This equation can be turned into a quasilinear one by taking a derivative with 
respect to <f> and cancelling the factors, since for the analysis of perihelion 
precession circular orbits are of no interest, 

«W = 5 *"(«)• (39) 

In order to solve these equations perturbatively we need expand them in 
powers of a small parameter. A good candidate is the inverse of the angular 
momentum per unit of mass, I, since according to Kepler's 1-2-3 law, which is 
assumed to be a good approximation at a great distance from the source, it 
behaves as / ~ y/rnr, where m is the mass of the particle. It can be combined 
with the mass of the source mo to yield an acceptable dimensionless small pa- 
rameter for analysing the far gravitational field. Hence we shall use e = m /l 
and expand u in the following way, 

n 

u = e 2 ^u„e" + 0(e 14 ). (40) 

n=0 

The reason for starting the expansion at this order is that the expression for 
the Kepler ellipse, which is expected to be the first term, is second order in e. 
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The energy per unit of mass of the particle is also to be expanded in e. 



Therefore we write, 

n=n 

E = l + e 2 J2 E n e n +0(e u ). (41) 

n=Q 

In order to avoid the appearance of secular terms we use a coordinate ip 
related to by, 

ijj = u)(j) uj= ^Jl + ^uoit 1 . (42) 

The coefficients c n are all of the order e 2 except C2 which is clearly of zeroth 
order in e and therefore equation (|38|) takes the form of a hierarchy of forced 
harmonic oscillators which can be solved iteratively up to the order of accuracy 
provided by our knowledge of the metric, 

u n ^+u n = f n (ip). (43) 

The first terms of the expansion of the solution to the equations Ij38(l and 
(jSHl) are, 



= — (1 + v / l + 2£ cos^) (44) 
Too 

ui = (45) 



u 

to 



6 + 4£ 

u 2 = (46) 

too 

u 3 = (8 + 4E ) 1 -^. (47) 
m o 

Of course the term of lowest order is the Kepler ellipse if Eq is negative. 

From the information we have about the metric we can calculate u> up to the 
eleventh power of e. These terms are just what we need to calculate the expres- 
sion for the perihelion precession, so we shall focus on them. The expressions 
of the terms u n are not needed and therefore we shall no enclose them here. 

Instead of writing the results as a function of the integration constants m.; , 
it will be more useful to write them in terms of the Geroch-Hanscn multipole 
moments, Pi, the physical interpretation of which is more appealing. Bear in 
mind that the odd multipole moments are imaginary and have to be multiplied 
by — i to obtain the usual real expressions J n . To calculate these moments we 
shall make use of the procedure described in 0] . 

If we have and expansion of the Ernst potential £ on the symmetry axis in 
terms of the Weyl coordinate z, viz. 

oo 

£(p = 0)^C„z-(" +1 ) (48) 

n=0 
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then the multipole moments can be calculated as follows, 

P n = C n n < 3 (49) 

Pi = C 4 + j C (C? - C 2 Co) (50) 

P 5 = C 5 + 1 Co (C 2 Ci - C 3 C ) + ±& {Cl - C 2 Co). (51) 

As a function of these multipole moments the first coefficients in the expan- 
sion of the energy per unit of mass, E, read, 

E x = (52) 
E 2 

E 2 = -6-lO.Bo-^ 1 (53) 

E 3 = -(8 + UEo) 1 -^ (54) 
"o 

E 4 = -^-20Eo-13E 2 + ^ + (2 + 3E Q )^ (55) 

E 5 = -(56 + 104^ + 56^)^, (56) 



where Eq, the Keplerian energy, is a free parameter which has —1/2 as a lower 
bound, corresponding to a circular orbit. 

The frequency u is different from one and therefore the orbit is not closed. 
Between two consecutives perihelion approaches the test particle traces an angle 
2tt/oj. Hence the perihelion has shifted an angle, A(f> given up to the eleventh 
power of e by the following expression, 

A(f> = 2 7 r(o;- 1 -l)=7r{Ao + Ai+A2 + A4 + A8+ 

+ A 16 + A 32 + A 2x4 + A 2><8 + A 2><16 + A4 ><8 }, (57) 

where the shift has been split into different terms according to their origin: The 
first one, Ao, comprises the Newtonian contribution to the precession, that is, 
the terms which remain after taking the classical limit c — > oo. Of course only 
the gravitational moments are present since rotation has no influence whatso- 
ever in Newtonian dynamics. Since the speed of light, c, and the gravitational 
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coupling constant, G, have been taken to be one, the terms look rather alike 
in magnitude. However, if the respective factors are written (a factor G/c for 
each e, a factor G/c 2 for each P 2ra , a factor G/c 3 for each P 2ra+ i, a factor c _1 
for each I and a factor G~ 2 for each E Q ), the actual magnitude of every term 
is recovered. For instance, the first Newtonian term has a factor G 2 and the 
second a G 4 . For oblate gravitational sources the quadrupole moment, P2 is 
negative, whence it contributes to a positive shift of the perihelion in the first 
order. In the next order the shift contribution of the quadrupole is, however, 
always positive. (N.b. Eq, though negative for bounded orbits, has a lower limit 
about —1/2 which does not allow it to overcome the energy-independent term.) 
On the other hand the sedecimpolc P4 term bears the opposite sign to the one 
of the quadrupole, 

3P 2 4 J/105 45£ \ P 4 A05 15£ \P 2 2 \ 8 

The second term, Ai, comprises the contribution to the perihelion precession 
due to a spherically symmetric mass distribution, i.e. the Schwarzschild effect. 
It can be calculated exactly in terms of elliptic functions and is of the order of 
G 2 /c 2 . The contribution of every order is always positive, 

Ai = 6 e 2 +(^ + 15£o) e 4 +(^ + 165£:o) e 6 + 



159105 16725 E Q 705 EA 
I H - 

32 8 8 



e 8 + 



/ 1701507 216375 E , 20115 Ef 
\ 32 + \ 



| zurio rjQ \ 10 



(59) 



In the term A 2 we have included the influence of a dipole of rotation on 
the perihelion shift. It is of the order of G 2 /c 2 . Since its lower terms in e are 
odd, it is sensitive to whether the probe rotates in the same direction as the 
source does or not. It is positive if the angular momentum of the source and the 
orbital angular momentum of the probe are antiparallel and negative otherwise. 
This has a qualitative explanation in the fact that the perihelion shift due to 
a mass monopole decreases with the angular momentum of the test particle; if 
the source happens to be rotating and its angular momentum is J, then I is 
replaced in the first order by / + 2 J/r as it can be seen in equation l|38|) after 
substitution of P\ by i J, whence the 'effective' I increases if both momenta are 
parallel and it would be expected that the perihelion advance diminishes. In 
contrast, the quadratic terms in Pi are independent of the direction of rotation 
and arc always positive whence they induce a perihelion advance. The cubic 
terms in P\ behave as the linear ones, 

A 2 = ^i^ + (168 + 48£o)^e 5 -(120 + 24£ )^4e 6 + 
Pq Po Po 
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+ (2562+ 1020 P + 36 Pq 2 ) 



iPi 7 



Pi 



€ — 



( 65607 44607 P 195 ET\ Pi 2 R 



V 16 ' 28 ' 4 J P 4 
+ ( (36046 + 17640 E + 1356 E% - 16 E$) ^ - (2048 + 672 Po) % ' 



/ 10256685 + 1320387 P + 118305 E$\ Pi ^ + 



?2\ v2 

V 112 28 28 J Pq* 

/ 3927489 569361 P 70659 Pg ^ A _ 3 A iPi 



V 8 + ^- i + ^^ + 174 ^ + 15 ^Jp7- 



/ 2735961 341339 P 27429 Pg \ i Pi 3 
V 28 + 7 + 7 J ~P^" ' ' 



li 



Under the name A4 the relativistic terms depending only on the quadrupolc 
moment, P 2 , and the mass are comprised. The first correction has a factor G 4 /c 2 
in front of it. As it was to be expected, it does not depend on the direction of 
rotation and, as its Newtonian counterpart, it is positive for oblate gravitational 
sources. The quadratic terms as a whole are always positive, 



A fqn + 49FM P * f 6 / 25383 28305 Po 375 P 2 \ P 2 8 

A 4 - -(90 + 42P )^3e - + —— + —— j — e + 



, / 2686203 + 503379 P + 80187 P^ 



V 112 ' 28 28 I P 3 





- l^ + ^Eo-^—^)-^}< w (61) 



The symbol A 8 stands for the corrections to the perihelion shift due to a 
rotation octupole moment P3. The first correction is of the order of G 4 /c 2 . They 
bear the same relation to the dipole terms as the quadrupolc to the monopole 
terms: The overall sign changes, 

A 8 = -{(30 + 24P )e 7 + (801 + 825 P + 138 P5 2 ) e 9 + 

+ (^1 + 16475 P + 4278 P 2 + 102 P 3 ) e 11 (62) 

The influence of the sedecimpole gravitational moment, P4, is included in 
Ai6 and does not show up until the tenth order of the small parameter. In 
non-geometrized units it is proportional to G 6 /c 2 , 



*16 



(- + 570£o + ^)_^, m 



The last multipole moment to be considered is the rotational trigintaduopolc 
moment, P5, and it is comprised in A 32 . It is of the eleventh order in e, 
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^32 



"» + 210 « + i^S)« e ". (64, 



Now we review the couplings among the different multipole moments other 
than mass. Up to the order considered there is no coupling between the gravita- 
tional moments higher than the mass (except for self-couplings), but there are 
rotation-rotation couplings and gravitation-rotation couplings. The first one to 
appear is the dipole-quadrupolc coupling, A 2X 4- It has a factor of G 4 /c 2 in the 
lowest order. If both angular momenta are antiparallcl and the source is oblate 
(P 2 < 0), then the contribution of the bilinear terms is positive. The quadratic 
terms in the quadrupole Pi are again positive if J and I are antiparallel. Finally 
the quadratic terms in the dipole are positive if the body is oblate, 



i 2 x4 = - (90 + 24 Eo) lPl f 2 e 7 - (3939 + 2127 E + 126 £ 2 ) lPl f 2 e 9 + 

Po Po 

2n r /Hm \ ip x p^ 



P P ( ( 1419 \ 

+ (2280 + 900 E ) -j^e w + | ( — + 303 E a \ 



P 







/■ 2825301 + 507992^ + 75765^ + A iP^ ■ n ^ 



V 28 7 7 J P 

The only rotation-rotation coupling up to this order is between the dipole 
and octupolc moments. It comes under the name of A2X8 & n d it is at least of 
the order G 6 /c 4 . Since it is linear on both rotational moments this term is not 
sensitive to the direction of rotation. It is of the tenth order in e, 



A 2x8 = (1068 + 972 E a + 120 E%) ^-Je 10 . (66) 

-Po 

The higher rotation-gravitation couplings involve the rotational dipole and 
the gravitational sedecimpole, A 2x i6, and the rotational octupole and the grav- 
itational quadrupole, A^g. Both appear first in the eleventh order of pertur- 
bation and are at least of the order of e 11 . The term A 2x i6 is positive if P4 is 
positive and the angular momenta are antiparallel, 

/4005 ^ 135 i£\ iPiP 4 n 

A 2x i6= + 495£ + ^-^ J -j^Tt ( 6? ) 



A 4x8 = (^ + 348^ + 45^) ^e 11 . (68) 

It would be of great interest to know the range of applicability of this per- 
turbative expansion. In an appendix at the end of this paper a simpler case is 
studied: It is shown there for which values of the parameters the expansions are 
acceptable for the Schwarzschild metric. 
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IV Precession of the line of nodes of a nonequa- 
torial orbit 

Let us consider now a bounded orbit slightly departing from the equatorial plane. 
If the mass distribution of the source were spherically symmetric, the equatorial 
plane would not be at all privileged and the orbit would always intersect it 
at the same nodes. However this is no longer the situation when the source 
is not exactly spherical, since then the nodes precess due to the perturbation 
generated by higher order multipoles. In classical nonrelativistic mechanics the 
first contribution to the precession of the nodes arises from the quadrupole 
moment, as it is shown in an appendix at the end of this paper, whereas in 
general relativity it is the rotational dipole moment the first one to contribute. 

In the approximation used in this section we consider a geodesic on the 
equatorial plane and calculate the evolution of small deviations from it. Since 
the geodesic equation reads, 

x p + r£ CT x p if = 0, (69) 

it is straightforward that nearby geodesies which deviate from the original 
geodesic by a vector 5^ fulfill, 

5" + 2T p a d p x <J + T^ v x p x a 5 l/ = 0. (70) 

Since we are interested in small deviations from the equatorial plane, we shall 
focus our attention on the 9 coordinate. Taking into account that the reference 
geodesic lies on the symmetry plane 9 = it/ 2 and that the first derivatives of 
the metric with respect to 9 vanish on it, the geodesic deviation equation for 5 e 
reduces to, 

S e -^g ee g pa ,eex p x°S e = 0. (71) 

Instead of considering an arbitrary bounded reference geodesic on the equa- 
torial plane, we shall restrict ourselves to geodesic circles. Of course many 
interesting features will be lost, but we deem this paper would become even 
longer if we take them into account. 

In order to compute the evolution of the nodes with respect to the azimuthal 
angle on the geodesic, the previous equation needs to be divided by <j) 2 , which 
is constant on the circle. From now on we shall write 5 instead of 5 e to avoid 
cumbersome notations, 

6^ + n 2 5 = Q n 2 = - -i- g m g paM x p x\ (72) 

2 q) z 

where every function is to be calculated on 9 = tt /2 and r = R, the radius of 
the geodesic circle. 

The previous expression states that the nodes of nearby geodesies are sepa- 
rated by regular intervals of the coordinate <j>. If is different from one, these 
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nodes will travel around the geodesic circle instead of remaining at constant 
values of the azimuthal angle. Since it would be of interest to write down the 
result in a coordinate independent expression, it will be required to make use of 
the geodesic equations to remove the dependence on the radius of the circle, R, 
and also to cast the energy, E, as a function of the angular momentum of the 
circular orbit, I, so that the final result be a function of the multipole moments 
and I only. 

As it is well known, the timelike geodesic equations can be obtained from 
the Lagrange equations applied to a functional L — g^ v x v and the constraint 
L = — 1. Two of these equations, namely (|34|) and ((35(1 . have already been used 
to remove the dependence on the derivatives of the ignorable coordinates. The 
equation corresponding to variations of 9 is automatically satisfied. We are left 
just with two equations, corresponding to the constraint and the variations of 



l = f"/^ P3) 

p r} . 2Ef '^A drA . f (l^l dr (^) = 0, ,74) 

where we have used r = and 9 = tt/2. 

From the last equation we get E/l as a function of R, 



E _ -b±v / b 2 - 4ac 
T ~ ~2~a ' 



(75) 



= A 2 f 2 i?- 4 d r (r 2 f- 1 ) - r 2 d r f-2Af R- 2 d r A = 

= -2 P a R-' 2 + 0(iT 3 ) (76) 



b = 2Rr 2 fd r A-2Af 2 R- A d r (r 2 f- 1 ) = 

= -12iP 1 R- i + 0{R- 5 ) (77) 



c = f 2 R- 4 d r (r 2 f- 1 ) = 2R- 3 + 0(R- 4 ). (78) 

Since the energy must be positive, the solution with the minus (plus) signus 
in front of the square root equation corresponds to a positive (negative) /. On 
the other hand equation 1)73(1 furnishes I in terms of E/l which has obviously 
the same coefficients which were obtained for E in the previous section after 
substitution of Eq = — 1/2. 

The frequency of the precession of the nodes can be written now indepen- 
dently of the choice of coordinates. As was to be expected, the first term is equal 
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to one, corresponding to the frequency of the nodes when the mass distribution 
of the gravitational source is perfectly spherical. 

The contributions to the precession of the line of nodes of a timclike geodesic 
which is slightly tilted with respect to a geodesic equatorial circle have been 
classified in the same way as it was done for the perihelion shift: A classical term 
plus the relativistic terms, divided according to their multipole content. In order 
to write down the expressions in non-geometrized units, factors including G and 
c have to be included as it was done in the previous section. The description 
which was made there of the necessary factors for each term in A</> for the 
perihelion shift also applies. The information that we have of the metric allows 
us to calculate the coefficients up to Z~ 13 . 

There is, of course, no contribution from a pure mass monopole. Were the 
frequency equal to one, then the nodes would remain at constant <fr. Therefore, 
A(f> — 2 7r (Sl — 1 — 1) describes the angle through which the line of nodes has 
precessed in one revolution of the reference circle, 



For oblate objects (P 2 < 0), the influence of the mass quadrupole amounts 
to a delay in the precession of the line of nodes with respect to the <j> coordinate 
on the circle of reference. Therefore the line of nodes does not precess in the 
same direction as the perihelion. This fact should not be confused with the 
precession of the angular momentum vector in time, which is of course positive. 
The contribution from the sedecimpole term is negative for positive P4 and 
the one from the sexagintaduopole is positive for positive P§. There is also a 
classical coupling between P4 and P2. The nonlinear terms in the quadrupole 
moment bear a positive sign, 



Most of the dipole-dependent terms are sensitive to the direction of motion 
of the test particle relative to the rotation of the source. If the angular momenta 
are parallel (I has the same sign as J = — i Pi), then the linear and cubic terms 
in Pi induce an advance of the line of nodes. On the other hand the quadratic 
and quartic terms in the dipole moment are not sensitive to that relative sign 
and their influence always amounts to a delay, 



A<j) = 7T {A + A 2 + A 4 + A 8 + Ai 6 + A 32 + A 2 x4+ 

+ A 2X 8 + A 2X 16 + A4 X8 } 



(79) 




(80) 



A 2 = - 



4zP Pi 



18 iP Q 3 P 1 
T 5 



+ 



18P 2 Pi 2 



243iP 5 Pi 
2f 



+ 



4131P 4 Pi 2 



+ 
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+ f i96«W - 3861 ^) r» + 27294 f° 6p i 2 r" + 



4 7 7 
^ 37983 zP 5 Pi 3 _ 268515 ^P 9 Pi \ n 



+ I 7 32 



4o 4 , 1060043 P 8 ^i 2 \ ; -i2 



-2565P 4 Pi 4 + 22 , - j I 
( 734200 iP Q 7 P 1 3 4944807 zP n Pi \ . 13 

+ iv — ? 64 — y ■ (8i) 

The relativistic contribution of the linear terms in the quadrupole moment 
Pi is again negative for oblate sources. The quadratic terms furnish a negative 
contribution to the precession shift regardless of whether the source is oblate or 
prolate, whereas the classical quadratic correction is positive, as it was shown 
previously, 

24P 3 P 2 2799P 5 P 2 / 4 2 12396 P 7 P 2 \ 10 
A 4 - — —+ up +[-24P P 2 + jl + 

( 327223 P 6 P 2 2 361993 P 9 P 2 \ ,_ 12 
+ 308 + 2^ J Z ' (82) 

As it happened with the perihelion precession, the rotational octupole term 
has the opposite sign to the one of the linear dipole term and of course it is 
dependent on the direction in which the probe orbits, 



a 12 iP 3 P 3 , 156*P 5 P 3 , 3387 iP 7 P 3 17707 iP 9 P 3 

As - — v — + — p — + ^7^ — + ' (83) 

The relativistic sedecimpole term bears the same sign as its classical coun- 
terpart, although it is much smaller in magnitude. There is no relativistic term 
in P 6 , just the classical term which has already been described, 



120P 5 P 4 2905P 7 P 4 
Al6 = Jw 2712 ' ( 84 ) 

and again we have another permutation of sign; the term in P5 bears the same 
sign as the linear dipole term, 



45zP 5 P 5 1905 ? Po 7 P5 
A32 = 2^ ii*— (85) 

Now we discuss the coupling terms between multipole moments other than 
the mass. The bilinear coupling between the dipole and the quadrupole moment 
is positive for oblate gravitational sources which rotate in the same direction 
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as the test particle. The quadratic term in Pi is negative for oblate objects 
no matter in which direction they rotate. The quadratic term in P 2 is positive 
when the angular momenta J and I are parallel, 



. 18zP 2 PiP 2 44HP 4 P!P 2 273P 3 Pi 2 P 2 

A 2X 4 - Jj + ^9 JYo + 

/ 8HP 3 PiP 2 2 198369 iPo 6 PiP 2 \ n 3007003 P 5 Pi 2 P 2 
+ ^ 2 + 28 ) 308li2 + 

+ ( - 10899 ^ 5 ^ 2 _ US2iP ^P 2 + 550958^P^ ^ f _ ^ 

The rotational bilinear coupling between the dipole and the octupole moment 
is independent of the direction of rotation and it is positive when J = — i Pi 
and J 3 = — i P3 have the same sign. There is however a higher coupling which 
is quadratic in the dipole moment and contributes in the same way as A 8 does, 



* 186P 4 PiP 3 100201 P 6 Pi P 3 3168*P 5 Pi 2 P 3 
A 2X8 - ^ ^ ^ . (87) 

The last terms to be considered up to this order are the couplings between 
Pi and P4 and between P 2 and P3, which are both sensitive to the relative 
directions of rotation of the source and the probe particle, 



^2x16 



255 iP 4 Pi P 4 14193 iP 6 Pi Pa , qq . 
21^ 47^ (88) 



2HP 4 P 2 P 3 197HP 6 P 2 P 3 

a 4x8 = a! 15 ' ^ ' 

V Conclusions 

In this paper we have displayed the approximate general asymptotically flat 
stationary axisymmetric metric for the vacuum spacetime surrounding a com- 
pact source possessing a symmetry plane orthogonal to the symmetry axis. The 
calculations have been carried out for the Ernst potential up to the term in 
r~ 7 in the pseudospherical radial coordinate. This has been useful to calculate 
relativistic corrections to the classical orbits around a compact mass distribu- 
tion. In particular the perihelion precession on the symmetry plane and the 
precession of the nodes of a slightly tilted circular orbit have been calculated. 

Concerning the perihelion shift, it has been shown that the contributions of 
each gravitational and rotational multipole moment follow a curious pattern of 
alternation of signs: The term in the mass monopole is always positive, whereas 
the linear quadrupole term is negative for positive P 2 (the quadratic term in 
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P 2 is always positive) and the sedecimpole term is again positive for positive 
P4. The sexagintaquattuorpole Pq is outside the limits of our perturbative 
expansion. If we think for instance of an oblate source which is very close 
to a classical axisymmetric homogeneous ellipsoid, then we would also have 
alternating signs in the gravitational multipole expansion (P2 „ = 3 Pq (c 2 — 
<z 2 )™/[(2n + l)(2n+ 3)]), c and a being the lengths of the ellipsoid's semiaxis 
parallel and orthogonal to the symmetry axis, respectively) and every term 
would have a positive contribution. The linear and cubic rotational dipole term 
is always positive for counterrotating configurations of the source and the probe 
(the quadratic contribution is always positive regardless of the relative rotation) 
whereas the octupolc term is positive if J3 = — i P3 and I bear the same sign. 
The trigintaduopole term has the same sign as the dipole term. It should be 
pointed out that the energy dependent terms are not strong enough to affect 
these signs. The coupling between different multipoles preserves the sign of the 
product of the corresponding linear terms, that is, if there is a P n — P m coupling, 
then its sign is obtained from the product of the —1 factors multiplying the 
corresponding linear terms in P n and P m . No coupling between mass multipole 
moments of order higher than the monopole appears to this order except for 
self-couplings. The sign independent terms, such as the mass monopole term 
and the quadratic dipole and quadratic quadrupole terms, always give rise to a 
perihelion advance. 

The influence of the different multipole moments on the precession of the 
line of nodes of a particle departing by a small amount from the equatorial 
plane is somewhat different from the behaviour which has been shown for the 
perihelion precession: The alternating pattern of the signs of the linear terms in 
the multipole moments other than the mass is preserved both for gravitational 
and rotational moments independently up to the considered perturbation order: 
If the Pm were all positive, then the contribution of each Pm would be opposite 
to the one of Pm+i and a similar reasoning is valid for the rotational terms. 
However the relativistic coupling terms bear the opposite sign to the one which 
would be expected from the product of the —1 factors before the corresponding 
linear terms: For instance, if the P\ and P2 linear terms are positive then 
the corresponding bilinear coupling is negative. On the other hand, this rule 
is not valid for the Newtonian self-couplings. Another difference arises from 
the fact that the gravitational sexagintaduopole moment Pq does influence the 
precession of the line of nodes in this order of perturbation, whereas it does 
not affect the perihelion precession. Also, there were no couplings between 
the gravitational moments (self-couplings and mass-couplings excluded) in the 
perihelion precession, but there is one (the classical P2— P4-coupling) in the node 
precession. It is curious that, while the Newtonian sign- independent terms are 
always positive, the relativistic ones (the quadratic dipole and the quadratic 
quadrupole self-couplings) are negative. 

Of course most of these corrections are meaningless for astronomical pur- 
poses in our solar system, but are very likely to be relevant for highly rel- 
ativistic astrophysical objects, suchs as pulsars and blackholes, where other 
post-Newtonian effects [1] have been shown to be present. 
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A Schwarschild spacetime 

In this appendix we shall study the range of applicability of the perturbation 
expansion for bounded orbits in Schwarschild's spacetime. For this metric the 
geodesic equation (|38(1 can be solved exactly in terms of elliptic functions. In- 
stead of writing it as a differential equation for u as a function of <fi, we are 
writing it as a differential equation for cj). In this section u is no longer the 
inverse of the pseudosphcrical radius but of the usual Boyer-Lindquist radius, 

/ 2 P E 2 — 1 \ ~ 1 ^ 2 

<i> u = (2P u 3 -u 2 + -^u+— —) =g{u)- 1 ' 2 . (90) 

We are considering that E < 1 and therefore g(u) has at least one zero. 
For bounded motion we need three zeros so that the orbit ranges between the 
apsidal points. If we call these zeros a > b > c > 0, then equation l|90l) can be 
integrated JT] in terms of the elliptic integral of the first kind, F(j,q), in the 
region b > u > c, 



(0 - M J P ° ia n C) = F(7, 9) = f da (1 - q 2 sin 2 a)^' 2 (91) 



n 



T = arcsin y__ i = ^— 02) 

An expression for u in terms of (f> is easily obtained taking into account that 
the elliptic sine is the sine of F(j, q), 




(b-c)sn 2 { X l^-^^-M\. (93) 



Since the real period of the elliptic sine is iK(q) = 4F(iv/2,q), then our 
u function is 2 K (g)-periodic. Therefore the exact perihelion precession of the 
orbit of a test particle around a spherical non-rotating compact object will be, 

2VSg( g ) 
VPo (a - c) 

the perturbation expansion of which in e coincides with the one of Ai in equation 
JHH), as it was to be expected, 




(2n- 1)!! 



,1 



2 



q 2n • (95) 



The limits of the range of applicability of the previous expansion are q = 1 
(a = b) and q = (b = c). For both values, there are two zeros of g{u) which 
coalesce and the bounded motion is no longer stable. Therefore we are led to 
study the range of parameters for which g(u) has double roots. This happens 
when u takes either of the values u± which are solutions of g'(u) = 0, 
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1 ± VI - I2e 2 
u ± = ^ , (96) 

and therefore the allowed region is the one enclosed between the curves g(u±) = 
in the E 2 — q 2 parameter plane. This yields critical values for the energy per 
unit of mass, E 2 = 8/9, and the perturbation parameter, e 2 = 1/12. The 
perturbative approach is no longer valid beyond this point in the parameter 
plane since there are no stable bounded orbits. It is remarkable that for e 2 > 
1/16 there is not only a lower limit for the energy of the bounded orbit but also 
an upper limit. 

B Classical precession of the line of nodes 

In this appendix we shall briefly derive the classical expression for the precession 
of the line of nodes. 

The Lagrangian for the motion of a particle in a gravitational field is, 

L = i r 2 + i r 2 9 2 + i r 2 sin 2 6 <j) 2 - V{r, 9) (97) 

The equations of motion and conserved quantities which can be obtained 
from this Lagrangian are: 

E = l - f 2 + X - r 2 9 2 + 1 r 2 sin 2 9 ft + V(r, 9) (99) 
I = r 2 sin 2 9<f> (100) 
r = r9 2 +r sin 2 9 <j> 2 - d r V(r,9) (101) 

r 2 9 + 2rr9 = r 2 sin 9 cos 9 </> 2 - d V(r, 9). (102) 

Truncating the Lcgcndrc expansion at the sexagintaquattuorpole multipolc 
moment, Pq, we get the following expressions for the energy per unit of mass, 
E, and the radius, R, of a circular orbit on the plane 9 = it/2 in the far field 
region, 

p _ Pq PqJjl 3 P 5 P 4 + 9 P 4 P 2 2 , 14 . . . 

E -~2T 2+ ^ 8^5 ) ( 103 ) 



!_Po 3P 2 P 2 15 P 4 P 4 + 36 PjP 2 14 
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In order to obtain the oscillations about the equatorial plane of a slightly 
tilted bounded trajectory with respect to the circular orbit, we introduce a small 
variation in the equation i|lU2|) . The result will be divided by 2 to yield the 
evolution of 59 as a function of d>, 



(S9)^ = -n 2 S9 n=dl + -p-V ee (R,Tr/2). (105) 

From this expression we get the frequency of the oscillations, f2, which can 
be written in terms of I and the multipole moments by inserting equation (|104fl 
into it, 



+ »(105P^ e + 45^P 2 P 4+ 81P^|) +o(rlS)i (106) 

which obviously coincides with the term Ao which was calculated using the full 
relativistic theory in section ITvl 
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